Degenerate Bose-Bose mixture in a three-dimensional optical lattice 
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We produce a heteronuclear quantum degenerate mixture of two bosonic species, 87 Rb and 41 K, 
in a three-dimensional optical lattice. On raising the lattice barriers, we observe the disapperance 
of the inference pattern of the heavier 87 Rb, shifting toward shallower lattice depths in the presence 
of a minor fraction of 41 K. This effect is sizable and requires only a marginal overlap between the 
two species. We compare our results with similar findings reported for Fermi-Bose mixtures and 
discuss the interpretation scenarios proposed to date, arguing that the explanation may be linked 
to the increased effective mass due to the interspecies interactions. 

PACS numbers: 03.75.-b, 05.30.Jp, 73.43.Nq 



Quantum degenerate gases are formidable systems to 
shine light on fundamental quantum phenomena occur- 
ring at extremely low temperatures, such as supercon- 
ductivity and superfluidity. In combination with opti- 
cal lattices and scattering resonances, degenerate gases 
give rise to strongly correlated systems, enriching even 
further the breadth of phenomena that can be directly 
probed. Indeed, the pioneering experiment on the su- 
perfluid to Mott-insulator transition [l[ has shown how 
physical models long studied in the field of condensed 
matter can be realized almost ideally. With two differ- 
ent atomic species, the wealth of quantum phases grows 
to a daunting complexity [2[, only marginally explored 
by experiments. Actually, experiments with heteronu- 
clear mixtures in three-dimensional (3D) optical lattice 
have been performed very recently only for Fermi-Bose 
systems [H, Q, while Fermi- Fermi and Bose-Bose mix- 
tures are yet uncharted territory. The importance of 
mixtures in optical lattices is hardly overstated: asso- 
ciation of dipolar molecules mapping of spin arrays 
0, schemes for quantum calculation and implemen- 
tation of disorder [8] represent only a few major research 
lines that potentially will greatly benefit from such sys- 
tems. In particular, Bose-Bose mixtures seem well suited 
for all these purposes, provided that collisional losses are 
adequately suppressed. 

This work reports the first realization of a degenerate 
Bose-Bose heteronuclear mixture in a 3D optical lattice. 
Exploiting the large mass difference, we investigate the 
regime where one species lies well in the superfluid do- 
main, while the other exhibits the disappearence of the 
interference pattern usually associated with the transi- 
tion from a superfluid to a Mott insulator. We focus only 
on the interference pattern and do not probe the shell 
structure characteristic of a Mott-insulator for trapped 
systems We observe that the presence of the 

superfluid enhances the loss of phase coherence of the 
second species occurring for increasing lattice strengths. 



Surprisingly enough, we find that the effect is sizable 
(see Fig. [1]), not only for small fractions of the minority 
species, but even for marginal spatial overlap between 
the two. 



The experimental setup has been described earlier [11 
LL2[. Here we briefly outline the production of a dou- 
ble Bose-Einstein condensate (BEC), first demonstrated 
in [13]. We capture 2 x 10 8 87 Rb atoms and a variable 
number of 41 K atoms in the stretched \F = 2^rriF = 2) 
states in a Ioffe millimetric trap (MMT) [l2|, Q3] and we 
start evaporation by lowering the MMT depth from 5 to 



3.5 mK, with final harmonic frequencies (ta 







2tt x (16.8,202,202) Hz [l5[. We then cool 87 Rb only 
by microwave evaporation, while 41 K is cooled sympa- 
thetically. A single sideband, 10 dB below the carrier, 
selectively removes 87 Rb atoms in the \F = 2,rriF = 1) 
state, which would quickly deplete the 41 K sample. BEC 
is achieved first for 41 K at 250 nK and then for 87 Rb: con- 
densates contain typically 4 (2) x 10 4 87 Rb( 41 K) atoms, 
with Thomas-Fermi radii of 2.4(2.2) /im and peak densi- 
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FIG. 1: (Color online): Interference pattern of 87 Rb for s — 
6, 11, 16 (left to right). Upper row: 87 Rb only; lower row: 
87 Rb mixed with 41 K (not shown). Each panel represents the 
column density over an area 175 /im x 135 /im. 
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ties of 1.4(0.9) x 10 14 cm- 3 . 

In the presence of 2 x 10 3 41 K atoms clearly condensed 
(inverted ellipticity after expansion), since no thermal 
cloud can be detected, our sensitivity places an upper 
bound of 50% on the thermal fraction of this species: 
this sets an upper bound on the temperature of 73 nK. 
Since the evaporation is stopped at the same level also 
in the absence of 41 K, the same temperature boundary 
applies also for 87 Rb alone. At 73 nK, the peak density 
of the 41 K thermal atoms is lower than 1/60 that of the 
condensate. For this reason, we will throughout neglect 
the 41 K thermal fraction in the following. 

The clouds are only partially overlapped: due to dif- 
ferential gravity sag, the 87 Rb condensate lies 3.2 /am 
below 41 K. Also, numerical integration of the 3D Gross- 
Pitaevskii equation (GPE) shows that both density dis- 
tributions are deformed and almost completely phase sep- 
arated, due to the strong K-Rb repulsion [16]. Indeed, 
the interspecies scattering length ax-Rb = 163an [17[ is 
much larger than both aRb = 99ao [18[ and ax = 65ao 
[l^ . with ao the Bohr radius. 

Once double BEC is achieved, we ramp three retrore- 
flected optical lattice beams at = 1064 nm, with 
frequencies differing by tens of megahertz, propagating 
along the x, y, and z axes and focusing at the center of 
the MMT with waists equal to (90,180,160) /im. The lat- 
tice strength is calibrated by means of Bragg oscillations 
and Raman-Nath diffraction [20[, yielding a systematic 
uncertainty of 10%, which must be added to the statistic 
uncertainties quoted hereafter. We choose the duration 
and time constant of the exponential turn-on profile, 50 
and 20 ms, respectively, by maximizing the visibility of 
the 87 Rb interference pattern. We wait for 5 ms with the 
optical lattice at full power and abruptly switch off the 
lattice beams (< 1 fis) and the MMT current (~ 100 /xs). 
We image both species after a typical time of flight of 
15 to 20 ms with resonant absorption: the interference 
peaks of 87 Rb are progressively smeared as the lattice 
strength is ramped at higher values. 

We measure the width of the central peak and the 
visibility v of the interference pattern, defined as [22[ 
v = (N p - N d )/(N P + Nd), where N p and N d denote, re- 
spectively, the summed weights of the first lateral peaks 
and of equivalent regions at the same distance from the 
central peak along the diagonals [see Fig. [5(e)]. At 
our lattice wavelength, the optical potential is the same 
within 10% for 87 Rb and 41 K. However, due to the their 
larger mass and scattering length, 87 Rb atoms tunnel less 
and repel each other more: the combined effect is that 
87 Rb loses phase coherence at lower lattice power than 
41 K. 

The main experimental result of this work is the obser- 
vation that the lattice strength at which the 87 Rb inter- 
ference pattern starts washing out is greatly shifted not 
only by a minor admixture of 41 K atoms, but also with 
a marginal spatial overlap. These findings, displayed in 
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FIG. 2: (Color online) Visibility and width of the central 
peak of the 87 Rb interference pattern, for partial (a,c) and 
zero (b,d) overlap with the 41 K condensate. In each panel we 
compare data and fit (dashed lines) with 41 K atoms (red) and 
without (black), i.e. with N K = 2(1) x 10 3 and 0. The insets 
show: how visibility of 87 Rb is extracted from images (e) and 
the interference pattern of 41 K visible for s=20 (f). 



Fig. m augment those of similar experiments carried out 
with the mutually attractive 40 K- 87 Rb Fermi-Bose mix- 
ture 0,H. 

To quantify the shift of the transition point, we plot 
the visibility and the width of the central peak versus 
the lattice strength in units of the 87 Rb recoil energy 
s [2l|. We compare data taken with 41 K and when 
41 K was not even loaded in the MMT (Fig. [2]). We 
fit the visibility with a phenomenological Fermi function 
v — uo/[l+exp(a(s— s c ))], which has the expected flat be- 
havior below the critical lattice strength s c and decreases 
exponentially for s ^> s c . From the fit, we find that, for 
^Rb = 3 x 10 4 , the critical s value is s c = 16.8 ± 0.4 for 
87 Rb only, while s c = 12.4 ±0.3 with N K = (2±1) x 10 3 . 

Using the formula rj = 0.696s- 01 exp(2.07^s)(a/ X L ) 
[221 ]. we relate s with the ratio of interaction to tunnel- 
ing matrix elements, rj = U/(6J), in terms of the lattice 
wavelength and the scattering length a. From the 
fit values, we derive both the critical value r] c and the 
exponents of the visibility decay for s ^> s c , v ~ r]~ v \ 
without 41 K we find rj c = 12.3±|J and v = 2.3t£;^, but 
with 41 K we have rj c = 3.9±iJ and v = 3.4t°' 8 (error 
bars are dominated by the calibration uncertainty of the 
lattice strength). 

The width of the central peak measures the inverse 



3 



correlation length: it signals the transition with a stark 
climb at s ~ 10 in the presence of 41 K and s ~ 14 in 
the absence thereof [see Figs. [2fc) and[2fd)]. After the 
transition, the width increase is steeper in the presence 
of 41 K even for null overlap (Fig. Eld))- this behavior 
is unexpected and so far unexplained. While the transi- 
tion points detected by the width and the visibility are 
different, the shift is the same, As c c± 4. 

Such a shift is surprising, given the little overlap be- 
tween the two species. Numerical integration of the 3D 
GPE with an s = 11 vertical lattice shows that the 
overlap is restricted to one lattice site out of 11 (see 
Fig. [3]). Generalizing to a 3D lattice, we expect that 
approximately only 10% of sites are simultaneously filled 
with both 87 Rb and 41 K atoms. In order to check that 
the shift is genuinely related to the interspecies interac- 
tion, we repeat the experiment with larger differential 
sag of the two samples. Once the double-species BEC 
is achieved, we relax the MMT harmonic frequencies to 
Co = 2tt x (9.2, 108, 108) Hz, thereby increasing the verti- 
cal separation to 11 /im, so that the overlap is totally neg- 
ligible even for any undetected 41 K thermal cloud (1/e 2 
radius = 7.7 /im). Figures [21(b) and^d) show that, al- 
though the visibility is slightly lower in the presence of 
41 K, the transition point is the same within our error 
bars: s c = 15.8(0.5). 

As a further test, we investigate the effect of the lattice 
turn-on ramps. As mentioned, for 87 Rb alone, the visibil- 
ity is maximum for a ramp of 50 ms. A priori, this ramp 
duration could be insufficient in the presence of 41 K im- 
purities. Therefore, we monitor the visibility versus the 
ramp duration. The turn-on profile follows essentially an 
exponential function with a time constant equal to 0.4 of 
the total ramp duration, smoothed at the end. We find 
that the visibility decreases for turn-on times longer than 
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FIG. 3: (Color online): density distribution of 87 Rb (red, left) 
and 41 K (blue, right) atoms, in the plane z — 0, calculated 
by numerical integration of the GPE: N^h = 3 x 10 4 , Nk = 
2 x 10 3 , vertical lattice strength s = 11. Shaded area (green) 
highlights the overlap, while the inset plot reports density 
profiles along the vertical line through x — 0. 



50 ms also in the presence of 41 K (see Fig. H]) and indeed 
the visibility difference is roughly constant at 15%. We 
conclude that, while in principle one could expect that 
in the presence of 41 K the lattice ramp should be slower, 
there is no evidence that this is indeed the case, at least 
as far as the loss of visibility is concerned. Moreover, 
we do not observe a significant change in the number of 
atoms with the duration of the ramp; hence three-body 
losses play a negligible role, consistently with the small 
overlap of the two clouds. 

As mentioned, our experimental findings extend those 
of Refs. [1, 0] to bosonic impurities and help to shine 
some light on the induced loss of coherence, whose origin 
has so far remained unclear, if not controversial. Giinter 
et al. report that the transition point, extracted from 
the visibility data with a fitting procedure equivalent to 
ours and expressed in terms of r] Cl shifts from 6.5 to 2.5 
for an 8% admixture. Both values are lower than ours. 
Ospelkaus et al. instead quantify the shift in terms of s c , 
as we do, and fit the visibility with the Fermi function. 
The reported shift As c = 3(1.5) for a 7 % fraction of 
impurities agrees with our result. With the same data, 
but measured from the width of the central peak, the 
transition point shifts from r\ c = 9 to 4.5 [3] and by As c = 

i(i) 1. 

Several different scenarios have been proposed to ex- 
plain the impurities- induced loss of coherence. We dis- 
cuss here the impact of our experiment on these scenarios. 

First, we notice that one would naively expect that 
strongly interacting impurities would locally increase the 
87 Rb atom filling factor and hence the chemical potential, 
by expelling or attracting them. However, increasing the 
chemical potential leads to an enhancement of the su- 
perfluid fraction, which is indeed the result predicted by 
Ref. [5], but opposite to what we observe. Ospelkaus et 
al. argue that fermionic atoms act like randomly local- 
ized scatterers for the boson order parameter, thereby 
splitting it into isolated superfluid domains unable to 
maintain a single coherent phase throughout the sample 
0| • The suggested phenomenon, reminiscent of Ander- 
son localization, depends crucially on the fermionic na- 
ture of the 40 K impurities, leading to localization. How- 
ever, bosonic 41 K atoms are expected to localize less, not 
more, than 87 Rb. Therefore, the interpretation based on 
Andersonlike localization does not hold in our case and 
might be unnecessary for the Fermi-Bose experiments. 

The loss of phase coherence could follow from heating 
of the two clouds when the lattice ramps up. In addi- 
tion to technical noise, which plays a role only at longer 
time scales, we consider the effect of the thermodynamic 
transformation involved. With two species, the effective 
masses change at different rates and thermal equilibrium 
requires an interspecies redistribution of entropy. In our 
experiment, this increases the entropy of 87 Rb and there- 
fore its thermal component. Quantitatively, due to the 
small overlap, we neglect the interspecies interaction en- 
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FIG. 4: (Color online) Visibility of the fringes at s = 12 versus 
the duration of the lattice ramp up. 



ergy and calculate the total entropy as the sum of two 
independent contributions. At s = 0, we calculate the 
contribution of each component from the formula for the 
total energy: E/N = 3/c B T c C(4)t 4 /C(3) +/i(l - 1 3 ) 2 / 5 (5 + 
16t 3 )/7, with t = T/T c , where T c and \i denote the crit- 
ical temperature and the chemical potential [23[. At 
s = 20 the contribution of 41 K is given by the above 
formula (replacing the atomic mass with the effective 
mass of the lowest energy band), while for the contribu- 
tion of 87 Rb we take that of a deep Mott insulator [24| : 
Sshe\i/k B = 327r 3 k B TY Jj Rj/(3rnu; 2 \l), where Rj and 
m indicate the radii of the superfluid layers of the Mott 
shell structure and the atomic mass. For an initial tem- 
perature of 73 nK ( 87 Rb thermal fraction equal to 0.21), 
the presence of 41 K increases the 87 Rb entropy and hence 
its thermal fraction, by 20%. As a consequence, the 87 Rb 
condensate fraction decreases from 79% to 75% and the 
visibility by approximately the same amount [25]. We 
conclude that this is a minor effect. 

A plausible interpretation of our findings invokes the 
dressing of 87 Rb atoms by 41 K: naively, dressed 87 Rb 
atoms weigh more, therefore tunnel less. This argument 
is made quantitative by Bruderer et al. [26] , who analyze 
the generation of polarons in a Bose-Bose mixture, i.e., 
quasiparticles composed of localized atoms ( 87 Rb) plus 
a cloud of superfluid phonons of the other species ( 41 K). 
The polarons have two different effects: at zero temper- 
ature they reduce the tunneling rate of 87 Rb atoms ac- 
cording to the picture outlined above; at finite tempera- 
ture, 41 K phonons undergo incoherent scattering, further 
hampering the hopping of 87 Rb. According to Ref. [261 ]. 
with the parameters of our experiment at s = 10, the 
energy scale relevant for the onset of the incoherent scat- 
tering is E p /kB ~ 1 nK and therefore it is comparable to 
and lower than the estimated lattice temperature, rang- 
ing from 10 to 30 nK, while the tunneling suppression 
at T = appears negligible. The analysis carried out in 



Ref. [26[ relies on suitably weak interspecies interactions. 
In our experiment the effective interspecies interactions 
are reduced by the low density of 41 K in the overlap re- 
gion but the approximations of Ref. [26[ are only weakly 
satisfied. The order of magnitude of the above value, 
however, proves that the inhibition of tunneling due to 
incoherent scattering of phonons deserves a deeper anal- 
ysis. 

Finally, we address the question of why the shift is sub- 
stantial even with the small overlap shown by Fig. 03 We 
conjecture that, owing to the inhomogeneous density, the 
loss of coherence starts from the borders of the 87 Rb con- 
densate where the filling factor is lower; the overlap with 
41 K, although small, occurs in this critical region. Alter- 
natively, one could call on collisions occurring during the 
expansion: since we do not observe any degradation of 
the 87 Rb Raman-Nath diffraction pattern caused by the 
presence of 41 K, time-of- flight collisions between the two 
species seem unable to significantly alter the 87 Rb mo- 
mentum distribution, at least in the superfluid regime. 

In summary, we have created a degenerate Bose-Bose 
mixture in a 3D optical lattice and studied the influence 
of a minor admixture of superfluid 41 K on the loss of 
87 Rb phase coherence. We find that this loss is furthered 
in the presence of 41 K, even if the two clouds overlap 
only at their boundaries. Our results agree qualitatively 
with similar observations carried out with a Fermi-Bose 
mixture, notwithstanding the different statistics of im- 
purities and the opposite sign of the strong interspecies 
interactions. Following a recent theoretical analysis, we 
find that the activation energy of polarons is of the same 
order as the estimated temperature in our optical lattice. 
Our degenerate Bose-Bose mixture will also be studied 
in the double Mott-insulator regime, promising for the 
formation of stable molecules. To this end, Feshbach res- 
onances predicted at moderate magnetic fields are instru- 
mental to control the interspecies interactions. 
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